Flow fields in soap films: relating surface viscosity and film thickness 
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We follow the diffusive motion of colloidal particles in soap films with varying h/d, where h is 
the thickness of the film and d the diameter of the particles. The hydrodynamics of these films are 
determined by looking at the correlated motion of pairs of particles as a function of separation R. 
The Trapeznikov approximation [A. A. Trapeznikov, PICSA (1957)] is used to model soap films as 
an effective interface in contact with bulk air phases, that behaves as a 2D fluid. The flow fields 
determined from correlated particle motions show excellent agreement with what is expected for the 
theory of 2D fluids for all our films where 0.6 < h/d < 14.3, with the surface viscosity matching 
that predicted by Trapeznikov. However, for thicker films with h/d > 7 ± 3, single particle motion 
is faster than expected. Additionally, while the flow fields still match those expected for 2D fluids, 
the parameters of these flow fields change markedly for thick films. Our results indicate a transition 
from 2D to 3D fluid- like behavior occurs at this value of h/d. 



PACS numbers: 47.57.Bc, 68.15.+e, 87.16.D-, 87.85.gf 



I. INTRODUCTION 

The motion of a particle in a viscous fluid causes a flow 
field to be created, that is, fluid mass is displaced in a 
very specific manner around the particle. The nature of 
this field depends sensitively on the geometry, boundary 
conditions and dimensionality of the system in question. 
For instance, flow fields in 3D decay as 1/R in an un- 
bounded fluid, where R is the distance from the localized 
perturbation. On the other hand, in a 2D fluid such as a 
thin film, flow fields decay logarithmically with distance 
[E EL 9 ■ One example of a 2D system exhibiting such 
long-range behavior under certain specific circumstances 
is that of soap films 0, [H, Q . A soap film in its simplest 
form consists of a thin fluid layer of thickness h buffered 
from air phases above and below it by surfactant layers. 
The fluid layer has a 3D viscosity T/buikj and the surfac- 
tant layers a 2D surface viscosity ?7i n t @, However, 
because of the finite thickness of the fluid layer, and the 
2D nature of the surfactant layers, it is unclear whether 
a soap film should be regarded as a 2D or 3D fluid for 
arbitrary film thickness, bulk viscosity and surface vis- 
cosity. 

Two particle microrheology [To| provides a powerful 
technique to determine the flow fields in soap films. This 
technique has been used with great success to character- 
ize the rheological and flow properties of 3D systems such 
as biomaterials [ll[ , polymer solutions [H, and bio- 
logical cells [lU- To a lesser extent, it has also been used 
to determine the surface rheological properties of 2D flu- 
ids such as protein monolayers at an air water interface 
[lH . Briefly, this technique looks at the correlated mo- 
tions of particles in the system of interest. When a par- 
ticle undergoes motion in the system, either by thermal 
excitation [3, [H[ or by an external perturbation it 
creates a flow field in the surrounding medium. This flow 
field affects the motion of other particles in its vicinity. 
Therefore, by measuring the correlated motions of pairs 
of particles as a function of their separation R, the flow 



field in the system can be determined. 

In this manuscript, we look at the flow fields in soap 
films of different thickness h and bulk viscosity Tfouik, 
while keeping the surface viscosity r/ mt constant. We em- 
bed probe particles of size d in these soap films, and 
are able to vary the dimensionless parameter h/ d by 
over an order of magnitude, from h/d = 0.6 — 14.3 [16j |. 
In order to account for the effects of the surfactant- 
laden interfaces and the thin fluid layer, we use the 
Trapeznikov approximation [171 ] that models the soap 
film as an effective interface with an effective surface vis- 
cosity rjT — r)hulkh + 2?7i n t . The flow fields in these films 
are then mapped as a function of R, and compared to 
theoretical models of soap film flow. Excellent agreement 
is obtained between the theoretical models and the flow 
fields in all soap films, irrespective of the parameter h/d, 
using rjT- However, we find that single-particle motion is 
faster than expected for thicker films with h/d > 7 ± 3. 
This leads us to state that the hydrodynamics of the soap 
films transition from 2D to 3D-like behavior at this par- 
ticular value of h/d. 



II. EXPERIMENTS 

A. Preparation of soap films 

Soap films are prepared from mixtures of water, glyc- 
erol and surfactants that stabilize the interfaces. By 
changing the ratio of water and glycerol, the viscosity 
?7buik of the soap solution can be controlled. The sur- 
factants used in this study are obtained from a com- 
mercially available dishwashing detergent brand Dawn. 
Known quantities of this detergent are added to the wa- 
ter/glycerol mixture to create the soap solutions. The 
chemical formulation of Dawn is proprietary and hence 
cannot be determined; however, for all the soap films in 
this study we use the same amount of dishwashing deter- 
gent (2% by weight) to ensure consistency. 



2 



Fluorescent polystyrene spheres (Molecular Probes, 
carboxylate modified, d — 210, 500 nm) are added to the 
soap solutions. The soap films are created by dipping a 
circular stainless steel frame of diameter f mm into the 
solutions and drawing it out gently. The frame is then 
enclosed in a chamber designed to maintain relative hu- 
midity and minimize convective drift. The particles are 
then imaged with fluorescence microscopy to determine 
their motion in real space and real time. 



B. Particle tracking by fluorescence microscopy 

The tracer particles are imaged in a fluorescent micro- 
scope at a frame rate of 30 Hz, with a 20 x objective (nu- 
merical aperture = 0.4, resolution = 465 nm/pixel) used 
for large particles (d = 500 nm), and a 40 x objective (nu- 
merical aperture — 0.55, resolution = 233 nm/pixel) for 
smaller (d = 210 nm) particles. For each sample, short 
movies of duration ~ 30 s are recorded with a CCD cam- 
era that has a 640 x 486 pixel resolution, with hundreds 
of particles lying within the field of view. The movies 
are later analyzed by par ticle tracking to obtain the po- 
sitions of the tracers [11]. From the particle positions, 
we determine their vector displacements by the relation 
Ar(i, r) = r(t + r) — r(t), where t is the absolute time 
and r is the lag time. Any global motion is subtracted 
from these vector displacements to minimize the effects 
of convective drift caused by the air phases that con- 
tact the soap film. These vector displacements are then 
used to determine the mean square displacement (MSD), 
(Ar 2 (r)}, where the average is performed over all par- 
ticles and all times t. Correlated motions of particles 
[t| [lj| are also determined by looking at the products of 
particle displacements, which we describe in Sec. IHII in 
greater detail. 



C. Determination of soap film thickness 

The thickness of the soap films in this study range from 
h ~ 300 nm - 3 /im, close to the wavelength of visible 
light. This make the spectroscopic technique of optical 
interference the most viable option for determining the 
thickness of these films. Immediately after taking each 
movie, the film is transferred to a spectrophotometer and 
its thickness h determined from the transmitted intensity 
[20| . We briefly describe the details of this technique: 
Two light rays of the same wavelength passing through a 
thin film will interfere with each other. This interference 
will be constructive or destructive, depending on whether 
one light ray has traveled an integer or half-integer mul- 
tiple of the wavelength with respect to the other ray. The 
light transmitted through the film will have a minimum 
(or the absorption will have a maximum) when 

(2ncos(9) • h = (m + 1/2)A (1) 



where n is its index of refraction of the film, A is the 
wavelength of light, 9 is the angle of incidence (typically, 
9 = 90°), and m is a non- negative integer. If X p k,i and 
A p fe,2 are two successive maxima (or minima) in the trans- 
mitted light then the thickness of the film can be easily 
determined by the relation (20j 
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FIG. 1: Absorption spectrum of a soap film (50:50 wa- 
ter/glycerol mixture, 2% Dawn by weight) as a function of 
normally incident wavelength of light. From the peaks in this 
spectrum and Eqn. [2j the thickness of the film is inferred to 
be h = 510 nm. Inset: The time dependence of h for two 
different soap films prepared with 60:40 water/glycerol ratio 
and 2% concentration of Dawn. 

Figure [1] demonstrates how the thickness of a soap film 
is determined in practice. It shows the absorption spec- 
trum for a soap film (50/50 water-glycerol mixture with 
2% Dawn) that has been placed in a UV-Vis spectropho- 
tometer (Agilent Technologies, Santa Clara, CA). The 
spectrum shows multiple maxima and minima because 
of constructive and destructive interference as light rays 
traverse through the soap film. By substituting the val- 
ues of the two successive maxima \ p k,i = 409 nm and 
Apfc,2 = 576 nm in Eqn. [5] and n = 1.4 for a 50:50 water 
glycerol mixture, the thickness h — 504 nm of the soap 
film can be estimated. We then average the value of h 
obtained by repeating this process for all the successive 
peaks observed in the spectrum, giving us h — 510 ± 10 
nm for this particular soap film. 

The inset to Fig. [1] shows a time series of film thickness 
for two soap films comprised of a 60:40 water-glycerol 
mixture with 2% Dawn as surfactant. It is evident from 
the figure that both films thin rapidly over an initial 
time scale of ~ 500 s, but rapidly equilibrate to a quasi- 
steady state over longer timescales where relatively small 
changes are seen in the thickness. Care is taken in our 
measurements to ensure that the particle trajectories are 
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recorded after this initial transient period. This has the 
added advantage that convective drift of the tracers is 
also substantially reduced beyond 500 s, resulting in more 
accurate measurements of particle motions. 



III. THEORY 

A soap film is considered "thin" when the thickness h 
of the film is comparable to the particle size d. In fact, 
it has been shown that thin films behave as a 2D fluid 
[SUB]- This assumption can be justified by modeling the 
soap film, with its two interfaces and a thin fluid layer, 
as a single effective interface in contact with two bulk air 
phases [17| (see Fig. [2]) . The effective interface then has 
an effective surface viscosity, tjt, which is given by 



(3) 



where ?7i n t is the surface viscosity of the two surfactant- 
laden interfaces. One immediate consequence of this ap- 
proximation is that a tracer particle in a soap film must 
diffuse as though embedded in an interface in contact 
with bulk air phases. According to Saffman [2lj], this 
diffusion follows the equation 
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where 77 a ; r is the viscosity of the air phase and 
(ho 7 7 is Euler's constant. It is expected that r/s = t\t 
and indeed this has been found to be true for thin films 
[l6[ ; the limits of this for thick films will be discussed in 
Sec. El 
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FIG. 2: Schematic of the Trapeznikov approximation, where 
the entire soap film is approximated as a single interface in 
contact with bulk air phases. Reproduced with permission 
from [11]. 

The other consequence of the Trapeznikov approxima- 
tion is that the hydrodynamics of a soap film must mimic 
that of a 2D fluid. To probe the hydrodynamics, we look 
at the correlated motions of particles embedded in the 
soap film. Similar to the treatment in Q, we measure 
correlated and anti-correlated coupled displacements of 
particles to determine the 'eigenmodes' of particle mo- 
tion in a 2D fluid. We give here a brief description of 
the method used by di Leonardo et al. [H . They actively 
perturbed pairs of particles by means of optical tweezers 
in a soap film. The displacements of the particles from 
their mean particle positions are then determined. These 
displacements are related to strength of the trap and the 



mobility of the particles. The mobility of each particle is 
affected by the presence of the other particle, and there- 
fore contains information about the nature of the flow 
fields in the soap film. The two dimensional Stokes equa- 
tion can then be solved and the resulting motions of the 
particles decomposed into eigenmodes of motion. There 
are four such eigenmodes in 2D, given by X x ± and X v ± 
where x, y represent motion parallel and perpendicular 
to the lines joining the centers of the particles and the 
± represent rigid motions and relative displacements re- 
spectively. These mobilities are given by the equation 
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A„± = b[l ± 
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(5) 



where R is the separation between the particles, L is a 
characteristic length scale and b has units of mobility 
(kg _1 s). This length scale L demands some explanation; 
flow fields in 2D are long-ranged and tend to diverge be- 
cause of the presence of the logarithmic term in Eqn. [5] 
This divergence is cut-off due to the presence of a length 
scale L, which has many possible origins. Some of these 
include finite size of the film, inertial effects, and viscous 
drag on the interfaces from the surrounding bulk fluid 
phases (air). In the subsequent sections, we will look at 
a series of soap films with different material parameters 
(??buik, h) and discuss in detail the origin of this length 
scale L. A detailed derivation of Eqn. O can also be found 
in reference [J]. It is clear from Eqn.[5]that the mobilities 
of the particles split around a mean mobility b, with rigid 
motions (+) being favored and relative displacements (— ) 
being opposed. The hydrodynamic interactions between 
the particles are also governed by the fluid layer of the 
film, shown by the appearance of r^ u \\Ji in the equation. 
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FIG. 3: Coupled motions of particles embedded in the soap 
films. There are 4 possible eigenmodes associated with 
this coupling, correlated and anti-correlated motion parallel 
(D x ±) and perpendicular (D y ±) to the lines joining the cen- 
ters of the particles. 

Our approach to determine the flow fields in soap films, 
while analogous to di Leonardo's approach, has some 
important differences. We measure the correlated ther- 
mal motions of pairs of particles embedded in the soap 
films. There are four such eigenmodes in 2D, represented 
by D x ± and D y ± which correspond to the longitudinal 
and transverse components of coupled motion, with the 
± representing correlated and anti-correlated motion re- 
spectively (refer Fig. [3]). The correlation functions are 
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given by [22|, \2c 



D x± (R,r) = (-{Arl(T)±Ari(T)f5(R~RV))^ 
D y± (R,r) = (l[ArUr)±Ari(r)} 2 6(R-R i %^ (6) 



where i, j are particle indices, the subscripts x and y 
represent motion parallel and perpendicular to the line 
joining the centers of particles and R lJ is the separation 
between particles i and j. The average is performed over 
all possible pairs of particles with a given separation R. 
This has the advantage of averaging out correlations from 
other particles that are not part of the pair. We are there- 
fore confident that many-body effects due to the presence 
of other particles are minimized in our correlation func- 
tions. Similar to [l5j . we observe that D x ±, D y ± ~ r 
which enables the estimation of four r- independent quan- 
tities {D x ±/t) t and {D v ±/t) t depending only on R and 
having units of a diffusion constant. In this manuscript, 
we shall discuss only these r-independent quantities. 

These r-independent correlation functions have a phys- 
ical interpretation; for instance, {D x+ /2t) t is the diffu- 
sion constant of the center of mass of the particle pairs 
along the line joining their centers, while {D x — /2t) t is 
the diffusion constant of the interparticle separation. Our 
correlation functions can then be trivially related to di 
Leonardo's eigenmobilities by the relation 



(D x± /t) = 2k B TX x± 
(D y ±/r) = 2k B TX y± 



(7) 



From Eqn. O we can then derive theoretical expres- 
sions for our thermally driven correlation functions which 
are given by 



(D x± /r)=B[l±CH-)] 



(D y± /r) = B[l ± C(ln( 



R' 



1)] 



(8) 



where B = 2k B Tb has units of a diffusion constant and 
C = l/(47r?7buik^fr) is a non-dimensional constant. In the 
subsequent section we explore the validity of these theo- 
retical expressions for a range of soap films with varying 
h/d and viscosity ?ybuik of the fluid layer comprising the 
soap films. 



IV. FLOW FIELDS IN SOAP FILMS 

Particle motions in the soap films are quantified by 
measurements of the mean square displacement (MSD), 
(Ar 2 ), which is ensemble-averaged over all particles in 
the field of view. Figure 0] shows the MSD for one par- 
ticular soap film (solid symbols, refer caption of figure 
for details) plotted against the lag time r. Also shown 
in the figure is the corresponding MSD for diffusion in 
a 3D solution that comprises the fluid layer of the soap 



film (open symbols). From the figure, it is clear that the 
MSD is linear with respect to r, indicating free diffusion. 
By comparing the two MSDs, it is also evident that dif- 
fusion in the soap film (solid circles) is faster than in the 
corresponding bulk solution (open circles). This makes 
sense, as the Trapeznikov approximation states that the 
particle is at an effective interface in contact with bulk 
air phases. Since the air phase has a significantly lower 
viscosity than the fluid layer, this speeds the diffusive 
motion of the particle. Finally, Eqn. 2] can be solved to 
estimate the effective surface viscosity rjs from the slope 
of the MSD. 
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FIG. 4: Mean square displacement (MSD) for a soap film, 
sample f in Table [I] (solid circles) compared to MSD in a bulk 
solution that comprises the fluid layer of the soap film (open 
circles). Dashed line shows a slope of 1, indicating free diffu- 
sion. The effective surface viscosity of the film can be evalu- 
ated from Eqn. [4] and is given by 7/ SjC ff = 8.03 nPa-s. 

After the effective surface viscosity has been deter- 
mined, we attempt to determine the flow fields in this 
soap film. This is done by looking at the correlated 
and anti-correlated motions of pairs of particles, as de- 
scribed in Sec. IIIBl of this manuscript. The four correla- 
tion functions (D x+ /t), (D x -/t), (D v+ /t) and (D v -/t) 
are shown in Fig. [5] as a function of particle separation 
R. From the figure, it is clear that the correlation func- 
tions split around a mean diffusion constant B, with rigid 
motions (+, solid symbols) being favored and relative 
motions (-, open symbols) opposed. Further, the cor- 
relation functions vary logarithmically as a function of 
R, evidenced by the linearity of the data on a log-lin 
plot. In fact, the data are well characterized by Eqn. [8j 
where B, C and L are fitting parameters. A visual 
way of determining the fit parameters is by noting that 
(D x+ /t) = (D x -/t) = B when R = L. Therefore, the 
dashed horizontal and vertical lines in the figure indicate 
that B — 0.54/im 2 /s and L = 64/^m, while the slope of 
the four correlations functions simply gives C = 0.13. 
These values have been used to fit the four correlation 
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FIG. 5: Correlation functions for the same soap film as in 
Fig- HI Symbols are: (D x+ /t), solid circles; (D x ~/t), open 
circles; (D v+ /t), solid squares; {D y -/r), open squares. Solid 
lines are fits to the data from Eqn |5J with B = 0.54/xm 2 /s 
(horizontal dashed line), C = 0.13 and L = 64/im (vertical 
dashed line). 



TABLE I: Material parameters for the nine soap films de- 
scribed in this paper. 7/buik (determined from diffusivity mea- 
surements in bulk solutions) has an error of ±5%, and values 
of h and d are certain to within ±2%. The uncertainties in 
7?int, derived from Eqn. 1 and 2, are given in the brackets. 



?7buik [mPa-s] 


h [rim] 


d [nm] 


h/d 


a. 2.3 


305 


500 


0.6 


b. 3.0 


640 


500 


1.3 


c. 6.0 


510 


500 


1.0 


d. 10.0 


1340 


500 


2.7 


e. 25.0 


1100 


500 


2.2 


f. 10.0 


780 


210 


3.7 


g. 25.0 


2184 


210 


10.4 


h. 30.0 


2100 


210 


10.0 


i. 30.0 


3000 


210 


14.3 
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FIG. 6: Scaled (a) longitudinal ({D x ±/t)/B) and (b) trans- 
verse {(D v ±/t) / B) correlation functions for five soap films 
plotted against the scaled separation R/d. Solid symbols rep- 
resent correlated motion while open symbols represent anti- 
correlated motion. Refer Table U for details about the soap 
films. Symbols are: triangles, sample a; hexagons, sam- 
ple c; squares, sample d; circles, sample f; stars, sample 
i. Solid lines are fits of the form 1 ± 0.13 ln(325d/i?) and 
1 ± 0.13(ln(325d/7?) - 1) respectively. 



functions with Eqn. [8] in Fig. [5] We should note that the 
length scale determined above is smaller than the field 
of view of the microscope, and is therefore not subject 
to finite size effects; that is, it truly represents a cut- 
off length scale for stresses in this particular soap film. 
Further, as B has units of a diffusion constant, it can 
be related to the self-diffusion of a single particle in the 
soap film (by replacing R = d in Eqn. [8|). Finally, the 
constant C — 0.13 represents how quickly the flow fields 
decay in this particular soap film. 

To test the validity of the functional form of the cor- 
relation functions, we measure and plot their values for 
different soap films with a range of bulk viscosities, thick- 
ness and tracer particle sizes. This is shown in Figs.^a) 
and (b), where we have non-dimensionalized the corre- 
lation functions by the scale factor B and the separa- 



tion R by the particle diameter d for five different soap 
films, including the one shown in Fig. [51 We find that for 
all soap films with h/d < 7, the longitudinal correlation 
functions scale onto a single curve, that is described by 
the equation 1 ± 0.13 \n(325d/R). For the film where 
h/d = 14.3 (stars), there is significant deviation from the 
scaling, particularly relating to the slope C (C — 0.06 
instead of 0.13). Similarly, the transverse correlation 
functions for the thin films can be well described by the 
form 1 ± 0.13 (ln(325d/i?) - 1), while the thicker films 
again deviate from the scaling (the transverse correlation 
functions are noisier than their longitudinal counterparts, 
hence the deviation is not as clearly visible) . The dashed 
horizontal lines in Fig. [Ha) and (b) depict the splitting 
of the normalized diffusion constants around a value of 1, 
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FIG. 7: Fit parameters for nine soap films, including the five 
shown in Fig. [6] as a function of h/d. (a) C; the dashed 
line indicates a constant value of 0.13 for all soap films with 
h/d < 7 ± 3. (b) L/d; the dashed line represents a value 
of 325, from the scaling of the correlation functions shown 
in Fig. [6] (c) B (solid circles) compared to the one-particle 
diffusion constant D s ,i p (open circles). The inset shows the 
ratio B/D St i p to be nearly constant with a value of 1.9 ± 0.1. 



while the vertical lines are placed at values of R/d = 325 
and 325/e (= 120) for the longitudinal and transverse 
correlation functions respectively. 

The fitting parameters from the five soap films de- 
scribed above, as well as four additional ones are shown 
in Fig. [71 The fit parameter C is plotted against h/d 
in Fig. EJa), where we see that C ks 0.13 for all films 
with h/d < 7 ± 3, while C « 0.06 for thicker films. On 
the other hand, the parameter L/d shown in Fig. |7J[b) is 
nearly constant for all soap films. Finally, in Fig. 0Jc) 
we plot both B and D s \ p as a function of h/d. Nei- 
ther the fit parameter nor the one-particle diffusion con- 
stant show any discernible trend with h/d; however, 
judging by their equidistant separation on a logarithmic 
scale, they seem to be correlated together in some fash- 
ion. Indeed, plotting B/D 3> i p against h/d, we find that 
B/D s ip w 1.9 ± 0.1 for all soap films. This is perhaps 
surprising, as D s \ p is a local measurement determined 
from the trajectory of single particles, while B is ob- 
tained from the correlated motion of pairs of particles. 

These parameters can be understood by considering 
various potentially relevant quantities which all have 
units of surface viscosity. The first quantity is T)b = 
7?t>urk^- This product shows up in the theory (Eqns. O 
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FIG. 8: The scaled fit parameter ksT '/ '(2ttBC) compared to 
the one-particle effective surface viscosity rjs (solid circles) 
and r\B (open circles). The dashed line indicates equality 
between the two quantities. 



[7J and |HJ), and is an appealing quantity because it de- 
pends on two physical parameters that are easy to mea- 
sure. The second quantity is the effective surface viscos- 
ity f)s as determined from the one-particle measurements 
(Eqn.[?]). The third quantity comes from the Trapeznikov 
approximation (Eqn. [3]) , which states that the entire soap 
film can be considered as an effective interface with a 
larger surface viscosity r] T = i] hu i k h + 2rj int = r) B + 2r/ int , 
thus suggesting that t]t is more relevant than the first 
quantity r/s- The fourth quantity is the combination 
k B T '/ '(2wBC), based on the two-particle correlations. 
From examining the theoretical expressions for the cor- 
relation functions described in Eqns. [5J and [SJ we 
see that B = 2ksTb and C — l/(47rr7buik^) ; so that 
UbT / (2nBC) = t\b- Thus, testing this equality is a test 
of that theory. In particular, the contribution to rjx from 
?7i n t has not been included in the expression for C, and 
so it is possible that the theory needs to be modified. We 
note that the fluid layer in the soap film studied by di 
Leonardo et al. [4] consisted primarily of glycerol, and 
presumably rji n t was irrelevant for their system, that is, 
r\B ~ Vt ^ Vint- This is not true for our soap films, as 
we have situations where t\b ~ ?7int as well as rjs 3> fyint- 
To test these conjectures, we compare these three 
quantities in Fig. [5J where we plot the scaled quantity 
ksT / (2itBC) against the measured one-particle effective 
surface viscosity rjs (solid circles) and t\b (open circles). 
For thin films with low surface viscosities (lower left cor- 
ner of the graph), we see the solid symbols for rjs are 
to the right of the open symbols for r/s, in agreement 
with the Trapeznikov approximation. From the differ- 
ence in these two quantities, one can extract the surface 
viscosity rji nt which we have done previously for these 
data, finding r]i nt = 0.97 ± 0.55 nPa-s-m [16]. However, 
for thicker films with higher surface viscosities (top right 
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corner of Fig. [8]), we see the opposite is true; the effec- 
tive surface viscosities rjs (solid circles) are lower than t\b 
(open circles), showing that Saffman's equation (Eqn.rjJ 
measures a viscosity that is too low. This occurs even 
though the thick films still behave as a 2D fluid accord- 
ing to their correlation functions, in other words, Eqns.[5l 
[7] and [8] still describe the two-particle correlations, albeit 
with different parameters. 

The comparison with ksT / (2ttBC) in Fig. [8] bolsters 
these conclusions. For thin films, 775 matches better 
with ksT/ (2irBC) (comparison of solid symbols with the 
dashed line, lower left of Fig. [5]). For thicker films, t\b 
matches better (comparison of open symbols with the 
dashed line, upper right of Fig.[8j. Our data suggest that 
in the theory leading to Eqns. 03 El and [51 we should re- 
place t]b = T]b u \kh with rjx] this minor correction would 
improve the results for thin, less viscous films, and be 
negligible for thicker viscous films such as those studied 
in Ref. |4|. The disagreement between ksT / (2ttBC) and 
•qs shows the breakdown of the applicability of Saffman's 
equation to extract a useful surface viscosity from one- 
particle data, for situations with h/d > 7, and in fact 
coincides with a breakdown of the equality rjs = tjt [IB] • 

To summarize, for all films we have fc^T / 1 (2-kBC) = rfr 
and for thin films we additionally have r\s =Tjr- The ap- 
proximation rjs ~ t)t is mathematically obvious for thick 
films, given the definitions of these two surface viscosi- 
ties. It is worth noting that our prior work suggests that 
it is necessary to measure rjs (in thin films) to be able 
to determine r\r, as the viscosity of the surfactant layers 
Tfrnt is usually not known ahead of time ; this current 
work suggests that measuring ksT / (2-kBC) is an addi- 
tional method to obtain r]T and thus r]- m t. 

Using these insights, we return to the data shown in 
Fig. [7] For low values of h/d, C is constant, as shown 
in panel (a); meanwhile, B drops, as shown in panel (c). 
This suggests that the mobility b scales as b ~ 77^ 1 for 
thin films. Indeed, this explains the correlations of B and 
D s i p seen in Fig. Etc) for thin films, as D St \ p ~ 1/j/t 
would be expected [16J. However, at thicker films, C 
decreases, suggesting that b grows faster than 77^ ; the 
mobility of particles increases. This is consistent with 
Saffman's equation incorrectly yielding a surface viscosity 
r)s that is too low for the thicker films. 

We also note that the parameter L/d is nearly con- 
stant for all soap films. This is very surprising, as the 
logarithmic cut-off in the decay of the correlation func- 
tions can arise from three possible factors: 1) finite size 
of the film; (2) inertial effects; and (3) viscous drag on 
the soap film from the bulk air phases. The size of the 
frame used to house the film is of order ~ I cm, which is 
too large to explain the values of L obtained from fits to 
the correlation functions. Inertial effects arise at length 
scales given by L = 7/buik/ pU, where p is the fluid density 
and U is a typical probe particle speed. This gives length 
scales of order ~ lm, which is again too large to explain 
our fit-derived length scales. Finally, the viscous cut-off 
length scale is given by L — 77 Sl eff/?7air- Clearly, as can be 



seen from Fig. [8j this length scale must increase as the 
effective surface viscosity increases. Therefore, the near- 
constant value of L j d for soap films over a range of 77 Sj eff 
remains a mystery. Note that none of these potential ori- 
gins for L would predict any dependence on the particle 
size d; that is, they might predict a film-independent L 
but not a constant L/d as we find (from data including 
two different particle sizes d differing by a factor of 2.4). 

One additional particle-based explanation for the con- 
stant value of L/d could be capillary effects. Capillary 
effects arise from deformation of the interface by the par- 
ticle, with an energy gain obtained when two particles 
stick to each other. However, it isn't obvious that this 
moderately short range interaction would lead to such a 
very long length scale as we observe (L/d ~ 325). We 
leave the reasons for the constant value of L/d as a mat- 
ter for future research. 



V. CONCLUSION 

We have used the technique of two-particle microrhe- 
ology to characterize the flow fields in soap films of vary- 
ing thickness h, where 0.6 < h/d < 14.3 (based on probe 
diameter d). In particular, we determine the 'eigenmo- 
bilities' of thermally correlated motions of probe parti- 
cles in these soap films. These eigenmobilities consist of 
correlated motion parallel and perpendicular to the lines 
joining the centers of pairs of particles, and rigid and rel- 
ative motion as well. The eigenmobilities are found to 
split around a mean value and decay logarithmically as a 
function of particle separation for all soap films. The flow 
fields of all films we observe are well described by theoret- 
ical models. A study of the fit parameters shows that thin 
films have a simple behavior, where a surface viscosity t]t 
predicted by Trapeznikov over 50 years ago correctly de- 
scribes both one-particle and two-particle motions [l7l ]. 
In fact, for the thinnest films, our results suggest that 
prior theoretical work should be modified slightly to use 
Vt 0- In our work we also a transition from 'pure-2D' 
to '3D-influenced' behavior at h/d = 7 ± 3. For thicker 
films, two-particle correlations still follow the predicted 
form for a quasi-2D fluid, with surface viscosity tjt- How- 
ever, one-particle motion is faster than expected, and 
other significant deviations from the thin-film behavior 
are noted. The results of our study on soap films can 
have important consequences for other 2D systems, in- 
cluding but not limited to microfluidic flow in confined 
geometries, protein and surfactant monolayers at an air 
water interface and the cell membrane. 
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